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Abstract 

We show how to reformulate Variable Speed of Light Theories (VSLT) 
in a covariant fashion as Variable Light-Cone Theories (VLCT) by intro- 
ducing two vierbein bundles each associated with a distinct metric. The 
basic gravitational action relates to one bundle while matter propagates 
relative to the other in a conventional way. The variability of the speed of 
light is represented by the variability of the matter light-cone relative to the 
gravitational light-cone. The two bundles are related locally by an element 
M, of 5-^(4, R) . The dynamics of the field M is that of a SL{4, i?)-sigma 
model gauged with respect to local (orthochronous) Lorentz transforma- 
tions on each of the bundles. Only the "massless" version of the model 
with a single new coupling, F, that has the same dimensions as Newton's 
constant Gn, is considered in this paper. When F vanishes the theory 
reduces to standard General Relativity. 

We verify that the modified Bianchi identities of the model are consis- 
tent with the standard conservation law for the matter energy-momentum 
tensor in its own background metric. 

The implications of the model for some simple applications are exam- 
ined. We discuss the Newtonian limit, the appropriate generalisations of 
the flat FRW universe and the spherically symmetric static solution. We 
conclude that the variability of the speed of light in the early universe is 
a possible homogenising mechanism. However an examination of the post- 
Newtonian approximation shows that on the basis of the results of VLBI 
measurements, the new coupling satisfies F/Gn < 3.2 x 10~^ . We point 
out that other "massive" versions of the model with different asymptotic 
properties may still permit consistency with General Relativity at large 
distances while predicting departures at short distances. 

DAMTP-99-102 
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1 Introduction 



Albrecht, Magueijo and Barrow have proposed Variable Speed of Light 

Theories (VSLT) to provide a stable basis for the observed isotropy and flatness 
of the universe. In these models the speed of light is parametrised as a kind 
of equation of state tied to the radius of the universe. They demonstrate that 
with the appropriate dependence of the velocity of light it is possible to create a 
cosmic dynamics containing long time attractors for the evolution of the universe 
that provide an explanation for its current state. However the formulation of the 
theories is not obviously covariant and the implication of the theories for regimes 
other than the early universe is not immediately clear. 

Moffat 0, 1^ earlier presented related ideas and Clayton and Moffat have 
proposed an ingenious bi-metric theory that overcomes some of the problems in 
formulating covariant VSLT. This is more fully dynamical than the bi-metric 
theory of Rosen 0, [ll|, IH . 

The possibility of "superluminal propagation" induced by quantum effects 
has been discussed previously |]13|-||231- The results are intriguing but remain 
controversial. 

In a spirit of exploration we propose in this paper a theory with variable speed 
of light that is both geometrical in structure and explicitly covariant in the same 
way as standard General Relativity. Furthermore the energy momentum tensor 
of matter is conserved in the conventional way. In this way it achieves some of 
the same aims as ref which appeared while this work was in preparation. The 
proposal is to introduce into the space-time manifold, two vierbein bundles. One 
is associated with matter and the other with gravity. The matter vierbein can 
be strained relative to the gravitational vierbein giving a geometrical meaning 
to variability of the speed of light. The propagation of light is of course associ- 
ated with the light-cone of the matter vierbein. Because of this unconventional 
approach to light and matter propagation we will refer to the Variable Light- 
Cone Theory (VLCT). It has a more general structure than just a varying speed 
of light. The theory in some aspects has the flavour of the bi-metric theories 

|T0|, 0, |12| but the emphasis is more on the relationship of the two vierbein 



bundles rather than on the two implied metrics. 

The dynamics of the theory is specified by giving gravity its standard cur- 
vature based action, Iq, and matter its standard action, Im, in the appropriate 
metric, and by introducing a linking action, I^, that controls the relationship 
between the two vierbein bundles. The full action, J, is the sum of all three 
terms, 

I = Ig + Il + Im . (1) 
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2 General Structure 



The theory is most clearly formulated in the vierbein formalism and makes fun- 
damental use of the resulting Lorentz gauge invariance. We introduce a vierbein 
bundle appropriate to gravity, {e^a}, with the associated metric 

QfJ-v ^fxa^v 1 (^) 

where raising and lowering of a-indices is carried out with the standard Lorentz 
metric rfab — {1, —1, —1, —1} ■ The inverse vierbein is {e"'*} so that 



e^„e- = 5; , e«%, = 5," , (3) 



and 

^A^^ = e'^^e^/ . (4) 

The vierbein associated with matter is {e^a} and the raising and lowering of a- 
indiccs is by means of the Lorentz metric, r^^g = {1, —1, —1, —1} . The associated 
metric is 

hrgf,u = e^se^ " ■ (5) 
The two vierbein bundles are related by a local hnear transformation 

e^a = e^„M% , (6) 

where we assume that the matrix M is an element oi SL{A,R) . This implies that 
the volume elements in the two bundles are the same. We denote the common 
value of the two determinants by J . 

J = det{e/} = dct{e/} . (7) 

It would be interesting to know if the assumption that det M = 1 could be 
relaxed. It seems to play a role in ensuring the conservation of the matter energy- 
momentum tensor. 

We denote the inverse matrix by M" ^ so that 

M\M\^5t , M\M^-,^5I . (8) 

We introduce vierbein connections for both bundles and associated coordinate 
connections. The relationship between them is achieved by requiring that the 
vierbeins are covariantly constant in the appropriate way. 

Df,e^a = d^e^a + ^^,a - rJj^eAa = , (9) 

and 

D^Cya ^ d^e^a + l^^ia ^^vb " "^lu^Xa = . (10) 
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The requirement that rjab and 77^5 be covariantly constant imphes that cu/^ab — 
—^ijba and U^g^ = -LO^la ■ 

It is convenient to define a covariant derivative of M that includes both the 
right and left vierbein connections, 

L'^M" , = d,M^ , + < , - M% d;/ , . (11) 

However in differentiating a second time the appropriate spatial connection, F^^, 
must be used. The covariant deriviative Z)^ can be extended in a similar way. Its 
effect on M is the same as that of but a second differentiation must use the 
spatial connection F^j, . 

Gravitational curvature tensors are defined so that 

rb 



and 



D« D«J \4 = i^.s^.y" , (13) 

where includes the left vierbein connection field, uj^ah but not the spatial 
connection, F^^ . Similarly includes the only the right vierbein connection, 
'^ixab ■ We have then 

Rabnv — d^LOi,ab — d^LO^b + <^^a ^fcb — ^ua ^l^cb , (14) 

with a similar definition for Rabf_iu ■ It follows that 

[D„ D^] M^.^R^ ,^,M' , - M« -,R' - 2Ci,D,M^ , , (15) 

where 

ci.-l{ri,-rt,) . (16) 

The quantity C^^ is the torsion tensor in the coordinate basis. 

3 Gravitational Action 

The gravitational action is the standard action 
where 

R^e'^'e'^'Rab^.u ■ (18) 

The vierbeins and the connections are treated as independent variables. If we 
vary the former then 

= ^ 1 d'xJde,, [e'^^R'' ^ - ^e^i?) . (19) 



The variation of the vierbein connection yields 
Switching to the fuU covariant derivative we get 
Using the result 

JD^V'' = d^iJV'') + 2JC^^V'' . 
we can integrate by parts and obtain finally 

If there is no other interaction in the theory we can deduce from the 
of these variations that C^^ = 0, and 

the standard equations for matterless gravity. 

4 Linking Action 

We arrive at the action for the matrix, M, that hnks the two vierbein bundles 
by treating it as the spin variable in a sigma model. This makes possible the 
gauging away of local (orthochronous) Lorcntz transformations in each bundle 
separately. The point is that only the local distortion of one bundle relative to 
the other is of physical significance and not the particular choices of local frame 
in the two bundles. 

Because it is gauge invariant under local (orthochronous) Lorentz transfor- 
mations in either bundle, the action that achieves this is 

h = / d'xJg^^Ti{j,j.) , (24) 

where F is a new gravitational constant with the same dimensions as G . The 
matrix valued current is given by 

j, = {D,M)M-' , (25) 

More explicitly 

j;'={D,M^-,)M'' . (26) 



(20) 

(21) 
(22) 

(23) 
vanishing 
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It is also convenient to define an alternative version of the current, appropriate 
to the barred vierbein bundle, j^**, as 



]f ={m-^D^mY =M\D^M''^ . (27) 

We treat the vierbein, which enters through g^^^ , the matrix, M, and the 
connections cj^at and oo^al independent variables. The result for the linking 
action from the vierbein variation is, 

5Il = / d'xJ5e,, [e'^^g'^^ - ^e^^''^) Tr(j^j,) , (28) 

Prom the left vierbein connection we have 

= ^ / d'xJdu^abf'"' , (29) 

and from the right vierbein 

5h = / d'xJ6u;^-,if'-- . (30) 

On varying the matrix M we obtain 

= / [mM-'){D,f - 2Ci^f)\ . (31) 

The quantity 5MM~^ is an arbitrary element of the SL{A, R) Lie algebra and 
is sufficiently general to identify the other factor in the trace which is also an 
element of the algebra. 



5 Matter Action 

We assume that matter is propagated in the vierbein background {e^g} ■ This 
seems a consistent approach since it implies that matter behaves in a conventional 
way in relation to the gravitational field it experiences. However the theory does 
change the relationship of this observed gravitational field to the distribution of 

matter density. 

If we assume that the matter action does not depend explicitly on the right 
connection oj^ali then the most general form of the variation of the matter action 
is 

SIm^- j d^xJ (Se^cT"" + Tr{SMM-^U)) . (32) 

However because we assume that matter behaves conventionally in the effective 
gravitational field we may write 

Wm = -^/ d'xJSg^^f^^'' , (33) 
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where T'^'^ — T"'^ is the symmetric energy momentum tensor for matter. Since 

Qixu = e^a^^ " , (34) 



it follows that 
However since 
it follows that 
and hence 



5Im^-J d^xJSe^ae, '^f'^" ■ (35) 
5e^, = 5e^„M% + e^„(5M% , (36) 
T"" ^ M^e^'^f"'' , (37) 



T\ = excT'^" = ex-ae, "T'^'^ = f'^'gxu • (38) 

If we adopt the convention that barred quantities, that is those approriate to the 
gravitational background of the matter, have spatial indices raised and lowered 
with the barred metric we can define 

T^ = f--^A. ■ (39) 
Hence we get the simple seeming result 

T\ = f\ . (40) 
However it is important to recall that 

In fact T'^'^ is not necessecarily symmetric. This does not cause any difficultry. 
We see also that 



(42) 



where [■ ■ ■] indicates the projection of the contained quantity onto the Lie algebra. 
Using the result ^ = e^^exa we can show that 



U^" ^e^^e/T''^-^r)^''T . (43) 



where T = T'' „ . Note that T = T^' = T . 



6 Equations of Motion 

We obtain the equations of motion by requiring that the variation of the total 
action is stationary. The result is 

S^i"''- - 8^^ (^0"^') - i^^Tta^,.)) - T", = . (44) 
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\- (D,f - 2Clf) +U = . (45) 



SttF 



JIM] = . (47) 



Eq(ll) implies that 



7 Bianchi Identity 

Just as in standard General Relativity it is necessary to check that the theory 
satisfies the integrability conditions associated with the Bianchi identity. In the 
presence of torsion these are changed to the following 

where ; /i indicates the covariant derivative . In the contracted version this is 
R^^-U^^r\ =R^\^Cp+2R\C^^^ . (50) 
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per ' p (J II 



If we set 



and 



g\ = r\-U;r , (51) 

F% = Tr (j-j, - , (52) 



p \j jp 2 p 
then eg (PI) can be expressed as 

We now take the covariant derivative, D^-, of this equation in order to test the 
relationship of the Bianchi identity to the appropriate divergence law for the 
matter energy-momentum tensor. Eq (|50|) gives immediately . From eg (|5^) 

we obtain 

F'^V = -.y) + ^"^'Tr (j'^(j.;, - 3,,^)) . (54) 
From eq (^5D we have 



Tr(j- ^^j^) = 2C^^^Tr(j-j^) - 87rFTr(f/j^) , (55) 
and using eq(^) we find 

Tiif .^^f) = 2^CPp^R'^' - IGTrFC^p^iT'^' - ^g'^'T) - SnFTiiUf) . (56) 



The cyclic property of traces allows us to write 

Tr(j^(j.;^ - J,;.)) = TrifilD,, D^]M)M-') . (57) 
If we introduce T^^o- and TZfj^a as Lie algebra matrices with the definitions 

{T^Hu)ab = Rahfia and {T^fia)ab = Rabfj.a 5 (58) 

then eq (|l5|) gives 

TrW^Ua;, - J,;.)) = MfH,^ - ^11,^ - 2Cp^^j^) . (59) 



Eg (^71) implies that the second term on the right yields zero under the trace. 
From eq(HU|) and eq(ES|) we obtain the result 



(60) 

Combining these results we find 

= -^^7'^'(2C;,/?%-C^',i?-^J-167rF(C'^^,r'^^-(7-^C;,r'^^)-87rFTr(f/j^) . 

(61) 

The covariant divergence of eq(E3) then yields 



T% = TriUj') + 2iC^^^T^'-g^^C^^,T'^^) , (62) 

alternatively 

T\.^^ = T:T{Uj,)+2iCP^T\-C;,T%) . (63) 
It is useful to note that 

3xab = -e,''{T\^-f\,)erb ■ (64) 

We have then 

Tr(f/j,) = -iV\^ - f I^)(e.,f/^'^e„ ^ . (65) 



Using eq(^3D we find 

Tr(f/j,) = -(r];^ - f gf ^ , (66) 

where we have used the result TrjA = . In turn this implies F^^ = f . We 
conclude that 

D^T\ = -{Tl^ - Tl^)f\ + 2{C;,T\ - C'^,T%) . (67) 
The relationship between the covariant derivatives and D^, implies 

D^T\ = D,T\ + (F^, - fl,)T\ - (F;, - f ;,)r% . (68) 



We obtain then 

D^T\ = {VI, - Vl,)T\ - (f;, - fL)T% . (69) 

However we have F^^ = f and find, replacing T"" ^ by T°";^, 

D„T\ = 2C:j\-2Cl^f% . (70) 

It is readily checked that this is equivalent to the standard conservation law 

D.f\ = Q , (71) 

where D^j is the covariant derivative incorporating the metric connection for the 
metric g^i, . We conclude therefore that the extended equations we have proposed 
for VLC gravity are consistent with conservation of the energy-momentum tensor 
for matter in the appropriate (barred) metric. 



8 Connection Structure 

One of the features of the theory is the role played by torsion. Here we show that 
the equations of motion do permit an analysis of the vierbein connections. 

From eq(|47|) we can obtain cu^^g in terms of the other variables. We can 
therefore eliminate uj^al from the equations of motion. We have 

= M„ %^m"^ , = %{,,5|M^ , . (72) 

We find then 

(73) 

This equation exhibits j^afe as a linear function of u^ab ■ Eq(il) relates the torsion 
tensor C^j, linearly to j^ab ■ The resulting equation determines uj^ab in terms of 
M"- ^ and e^a and their derivatives. 

From the covariant constancy of e^a we have 

T'^. = e^'{d,e.a + u;^Je,t) • (74) 

Hence ^ 

If we define 
and 

^^fiXu = "e^ ^ujf,ab , (77) 
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then we find 

^fiXu = Cxfj^i, — Cyfj_x + C^xy + uj^xu 1 (78) 
where Cj^xu is the metric version of oj^xv and is given by 

^^xu = ^ {e^ "df.exa - ""d^e^a + d^gx^, - OxQufj) ■ (79) 
Eq(^) can be expressed in the form 

Cafiv + C^^Qytj — C^j^g^fj = — Xo-^j, , (80) 



where 
and 

It follows that 



G 

Xafiu = -^jaM , (81) 
ja^u — 6^ 6^ jcrab ■ (^2) 

Finally we have solution 

^^^l\u = —Xxfiu + Xy^x — X^xu — X^xdfj.^ + X^^g^x + ^nXv ■ (84) 

9 Special Cases 

The above solution for uj^xv can be used rather straightforwardly in special cases. 
We consider a situation in which the vierbeins and the linking matrix have diag- 
onal forms, namely 

e^,a = Aar]^a and MV = A^5« . (85) 

It follows that 

ga/i^^-l^a/. M\ = A-l5f , (86) 

and that 

e^-a = KAari^-a and = A-^A-^rj^^' . (87) 

We enforce an obvious correspondence between the values of the symbols /x, a 
and a to give meaning to the t] and 5 symbols. In an appropriate sense we can 
write Aa = Aa = and = Aq = A^ where this is convenient. 



In this special case we can compute j^ab from eq([73|) to yield 



^Vab + -UJxab[l-^ 



b . 
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Hence 



and 



jA[a,6] = ^^AaJ 2 - - — 1 , (89) 



^A,. = ^5^5XA^Aa - ^ - ^) . (90) 



We also have 



= ^ [a-J - ^ [j-J ^'^"^^ ■ 

9.1 Flat Expanding Universe 

An application of the above equations is to a fiat expanding universe where we 
set 

Ao = l , A = A{t) , and Aq = A(t) , A, = A5(t) , (92) 

and i runs over orthogonal spatial directions. The metric is therefore 

ds^ = de - A\t)dv^ . (93) 

The barred metric appropriate to matter is 

df = A\t)dt^ -Al{t)A^{t)dr^ . (94) 

If we adopt the convention that a coordinate label n can take a time value 
which we denote by t or acquire a spatial character which we denote by the 
variables x, y or z, then rotational invariance tells us that the only non-vanishing 
components of ujx^i, are ujxty = —^xyt and that 

U^ty = ^8^y ■ (95) 

The same is true for the metric connection, and we find that 

ujxty = AA6xy . (96) 

The same considerations of rotational invariance imply that only X^ty = 
—Xxyt are nonvanishing and from eq (|90|) we find 

X - ^ u (2 (97) 



We find also 



X^l — -^^x'y'Xx'ty' ■ (9^ 
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From eq(p^ we have then 

G / 1 \ / \ • 

Uj5^y = — {^-2UjS:^y + 5^yA^ X ^1 U 5 ^1 yl ^ - " ] ^ ^^^Xy ■ (99) 

Finally 

A A 

= — , (100) 

where 

G ( A^ A|\ , , 

We will use this result later when we consider fiat FRW space-time. 

9.2 Static Rotationally Invariant Space-time 

Another application of the formalism is to the equivalent of the Schwarzschild 
solution. We will use isotropic coordinates, so that 

= A{r) , = B{r) , and Aq = A(r) , A, = Asir) . (102) 

The unbarred metric is then 

ds^ = A\r)dt^ - B\r)dr^ , (103) 

and the barred metric is 

ds^ = A\r)A''{r)df - A.l{r)B\r)dr^ . (104) 

Rotational invariance implies that only the components Uxty = —'^xyti ^ttx = 
—uJtxt and uJxyz may be non-zero. In the case of the metric connection we find 

OJxty = , (105) 

u;ux = AA'r, , (106) 
where "Tx is the x-component of the radial unit vector. We have also 

LJ^y;, = BE' {fy5^^ - fj^y) . (l07) 

The components of Xx^u that may be non-zero have the values, 

Xxty = J^^xty - ^ - , (108) 

X,.. = ^.,..(2-|-§) . (109) 
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and 

X,,, = . (110) 
Solving for cux/^i, we find a pattern similar to that for the metric connection, 

(^xty = , (111) 



and 
where 

and 



OJttx = -p-^o; , (112) 

G ( A2 k%\ , , 

^-^^ +4Fi^"A|" A^J ■ ^"^^ 

We will use these results in analysing the Schwarzschild-type metric in VLC 
theory. 



10 Weak Field Limit 

The weak field limit has the form 



S + V , (116) 



where 

Similarly we can set 
where 



e'^'le'^'l'' = V,. , (117) 
C/xa = + h^a , (118) 



e^e^'V^V,. ■ (119) 

The connections ujnab and tu^^g are first order quantities. The matrix M has the 
form 

M% = M(°)% + m% , (120) 
where M^^^ is a Lorentz transformation. We have then 

^(0) a^_^(0)^a _ (^21) 

Wc can use M and its inverse and e and e to convert superfixes and suffixes 
between the various bases. For example wc have 

m\ = m%M(°)% . (122) 
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The requirement that det M = 1 imphes that 

m\ = m% = m% = . (123) 

The relationship between e and e imphes that 

hf,a = hf,a + m^a , (124) 

together with corresponding equations in other bases. 
To lowest order 

-d^u^ab ■ (125) 

Hence 

= e(°)'^'^e(°)''^i?,fe^, , (126) 

so that 

= S^c., - - S.c.^ - , (127) 

and 

R = 2d^uj^'^'' . (128) 
Again in the lowest order approximation 

jfiab = df^rUab + UJf,ab " ^^/lab ■ (129) 

If we convert to the coordinate basis we have 

jfiXr = df,mxr + UJf^Xr " ^^fiXr ■ (130) 

We can also evaluate . In this lowest approximation it coincides with . From 
the equation of motion we have 

jtilx,T] = Jm[a,t] = . (131) 

That is 

dam[x,r] + ^^aXr " U^aXr = . (132) 

This result tells us that the torsion in the unbarred, gravitational vierbein bundle 
is zero. 

Eq(||) yields 

r]>'''d^{d^mxr + u;^Xr-i^aXr) = -87rFUxr , (133) 

where 

Uxr=fxr-\vXrf . (134) 

Making use of eq (|132|) we obtain the result 

V^''d^d^m^x,r} = -8nF (fxr - ^Vxrf'^ . (135) 
15 



The remaining equations are to lowest order 

- ^VaxR = SnGT^x , (136) 

and 

C',. = l{d,K'-duh^' + u^,\-^u\)=0 ■ (137) 

In the present approximation m[n^iy] can be absorbed by gauge transformations of 
the form 

LO^Xt ^aXr + 5<^0Ar , aud LO„Xt ^^aXr + da(j)XT ■ (138) 

We can assume therefore that in this approximation m^^^^] vanishes. Therefore 
nifj^^, may be assumed symmetric. It satisfies 

d^m^^ = -87iF{f^,~^ri^,f) . (139) 

The gauge invariance referred to above means also the we are free to choose /i^i/ 
to be symmetric with the result that h^u is also symmetric. 

Under these circumstances we can solve the vanishing torsion equation to give 

Wi^Am = d^Kx - dxK^, . (140) 

Eq(^^ now yields 

d.d^h^+d^d^h-' ^-d'h,,-d^d,h\-r],^{d,drh''^-d^h'' ^) = SnGT,, . (141) 
We now refine our coordinate system by choosing the harmonic gauge. 

9^X^ = . (142) 
In the lowest approximation it yields 

d,h\ = ^dxh\ . (143) 
The equation of motion then becomes 

- -9' f V - lv^..h^ r) = SttGT^. , (144) 



or 



d^h^^ = 8nG(f^^-^r]^^T^ . (145) 
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10.1 Static Case 



In a static situation 

92 = -V^ , and Too = T = p , (146) 
where p is the density of matter. We have then 

V^/ioo = 47rGp , and V^moo = 67rFp . (147) 

It follows that 

V%o=47rG^p , (148) 

where 

Gm = G+^-F . (149) 

We interpret /loo as the gravitational potential seen by matter. Therefore Gn can 
be identified with Newton's constant. Both gravitational constants G and F enter 
into the structure of weak gravity but only in a particular linear combination. 



10.2 Gravitational Waves 

A time-dependent matter distribution can act as a source of gravitational radi- 
ation. Clearly from eq( |145| ) we see that there are such waves of a conventional 
type. However eq (|139| ) shows that we may also have unconventional gravitational 
radiation associated with the oscillations of the field m^y . It is interesting to 
count the degrees of freedom introduced in this way. 

Dealing first with the conventional degrees of freedom represented by h^y 
the standard argument is that the conservation of the energy momentum tensor 
implies that the gauge condition eq( |143| ) is maintained. This is a necessary 
consistency check. Far from the source the wave, in the weak field limit, may be 
treated plane wave of the form 

V = Ve*'" > (150) 
and the gauge condition leads to the result 

k^A^ ,-hyA\ = Q . (151) 

These four conditions reduce the original ten degrees of freedom to six. Fi- 
nally four of these degrees of freedom may be removed by appropriate coordinate 
transformations that preserve the gauge condition. This leaves the standard two 
polarizations for the gravitational wave. 

When we come to the degrees of freedom represented by the m^y we see that 
the conservation law for the matter energy-momentum tensor implies 

d'^d^m\ = 2'KFd^T . (152) 
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However we know that 

f = ^O^h'.. (153) 

Therefore 

d^d^m\ = ^d''d,h\ . (154) 
This imphes, for waves that have their source in matter, that 

d^m\ - ^d,h\ = . (155) 

If we assume we are far from the source so that 

= B^^e^'^-^ , (156) 

then we find 



k^B\-—KA\ = . (157) 



F 

AG' 

The symmetric matrix m^j, has ten degrees of freedom. However one is removed 
by the vanishing trace condition. Four others are removed by the conditions 
above, leaving five remaining degrees of freedom. This imphes that there are 
seven degrees of freedom left in the combined {/i^,^, m^i,}-system. However the 
matter metric does not see one of these degrees of freedom. We can easily deduce 
the result that 

d^h\-l(l + ^)d,h\ = . (158) 



2G. 

These four equations reduce the original ten degrees of freedom in hf^^, to six. 
Because the coefficient of the second term departs from 1/2 when F 7^ 0, we 
cannot find a change of coordinates that removes further degrees of freedom. 
Matter interacting with the observed gravitational waves will see only six degrees 
of freedom. This is the conventional result for non-standard gravity theories. 



11 FRW Universe 

A strong motivation for the VLC theory is the possible modification that the 
theory may make to the evolution of the early universe. We examine this issue 
in this section. 

The equation of motion, eq (^5D , yields the result 

+ n^j' + ujf.f - fu;^ = -SttFU . (159) 
In the present case the only non- vanishing component of F^^ is 

rr, = ^ = 34 • (160) 
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In taking the {00} component of eq( |159| ) we find in the present case 



(161) 

The spatial components of the equation yield the same information. We find also 

and ^ 

T^^-'^*) = ^ fx) ' (163) 



yielding the equation of motion 



_i9.f4) Gif4V = 4.G(, + 3p) . (164) 



A \E} F3 V A 
The spatial components are 



'A\. I A^' 



i?% = -^9J-U.,-2 — 5., , (165) 



and 

Tr(r..)^-^(2-^-^j^^j 4. , (166) 



\ A2 A|\ / i ^ ' 



leading to the equation of motion 



These equations reduce to the standard equations for the fiat FRW universe when 
F = . They may be manipulated to yield the first order equation 

It is also possible to verify, consistently with the Bianchi identities, that they 
imply the the appropriate energy conservation law 

5*(^p)+p5t(^) =0 . (169) 

If we set 

c= — = A3 , (170) 
19 



then c is the velocity, seen from the gravitational background, of a signal traveUing 
with the velocity of light in the matter background. It is therefore the variable 
speed of light generated by the variable light-cone structure of the theory. We 
have 

The equations of motion can be rewritten in the form 
and 

11.1 Steady Expansion 

An interesting consequence of these equations is that if the pressure satisfies 

p = Ap , (174) 

where A is a constrant then there exist solutions with constant values of c 7^ 1 . 
We find 

( aV / , 1 



AE \ c 



and 



We have then 



- = 47rF(l + \)p , (175) 



For relativistic matter A = 1/3 and we obtain the quadratic equation 

'G 



(cf + 2(l-^)c^-3 = , (178) 



with the physical solution 



e^.e5^-(l-f)+2f-£ + g.l + | . (179) 

For pressureless matter A = and we obtain similarly 



3 / 2F\ 1 /36 / 2F\^ 44 3F , , 
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We take the relativistic case as an illustration. From eq( |176| ) we see that the 
expansion is controlled by 

(^)^5^,. (181) 

where Eq = E{cq) . Because the speed of light is constant we have from eq (|169|) 

Pocjj . (182) 

This is exactly like the standard scenario for expansion in the early universe with 
the change that Gn has been replaced by GEq . It follows that 

AocVt . (183) 

It is clearly true also that in the matter coordinates 

AocVi , (184) 

where A = Cq^^^A and t = c^^^t . This steady expansion is essentially the same 
as the conventional case. 



11.2 Non-steady Expansion 

It very difficult to investigate the solutions of eqs (|172D and (|173D in general. 
However some feeling for their nature can be obtained by examining solutions 
near the above steady expansion. If we set 

c = Coe« , (185) 

and assume that ^ is a small quantity and that we can neglect quantities quadratic 
in ^ and ^ so that eq( |173| ) becomes 

f A\^ SttG , , 

If we examine the case of relativistic matter, then when we combine this result 
with eq( |172D we get 

.A, / i \Vl / 2 l^ -^F 



Keeping only terms we obtain 
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where, in keeping with the approximation, A may be assumed to be the the 
steady expansion solution obtained previously. As a further simplification we 
will assume that F is very small and set F = in the above equation. We obtain 
then 

e+^e-|e = . (189) 
This homogeneous equation has two linearly independent solutions 

e = r , and ^ = Jf3 ■ (190) 



where a = (V33 - l)/4 = 1.1861 and (3 = (V33 + 1)/4=1.6861 . Using these 
solutions we can construct a range of scenarios for a (small) departure from the 
steady expansion. A solution of the form 

^ = ^ ' (191) 

will yield 

c = coe'/*' , (192) 

which for 6 > results in the speed of light increasing in the far past and tending 
to the steady expansion in the future. If 6 < then the speed of light will increase 
to the steady value from a smaller value in the past. If we choose a solution 

^ = ar (193) 

then when a > the speed of light will increase from the steady value or when 
6 < it will decrease from this value. By choosing arbitrary combinations of 
the two solutions any mixture of these scenarios can be obtained. Of course the 
analysis being perturbative, the results cannot be trusted beyond the point where 
^ ^ 1, so the ultimate behaviour at very small or large times remains unresolved. 



12 Static Rotationally Invariant Case 

We derive the equations for the static rotationally invariant system in isotropic 
coordinates. They can be used to study the equivalent of the Schwarzschild 
solution. However we will confine our attention here to the asymptotic properties 
of the solution as a way of investigating the post-Newtonian approximation. 
The equation of motion, eq (^5D , in the absence of matter yields 

d,f + T^,^f + uj,f-fcu, = . (194) 

Making use of the fact that 

Ia,- ^ - ^ +^ ^ ' ^^^^^ 
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we we see that 



and 



A 



) 

' B 



(196) 



Prom the {00} component of eq() we find 



A' 



'A'Y 2 A' B'\ A' _ 1 



(197) 

The other non-trivial components yield the same equation. Using the definition 
F^^ = Tr{j^>} we find 



" B\P \ Al A2 , 



and 



xy 



4 /tV^ 
HA. 

We also obtain for the relevant components of R/^i, the results 



Rtt 



A 
B 



'bp 



2 
r 



'bp 



- 2 



'bp 



Q_ 
B^ 



and 
where 



R 



xy 



■^xy "I" -^xy -^xy 



^xy 



B 



A' 

'bp 

Q_ 

S2 



+ 



Q_ 
B^ 



X' y : 



and 



xy 



A'\ Q 



IB 

r A 



A' 



rB^ 



Q_ 

S2 



^xy 1 



BV 



xy 



AP J B'^ rA\BP 
Prom the {tt} component of the equation of motion we obtain 



Vy) 



A' 



A' 



A' 



Q_ 

S2 



f[p 



2 - 



A2 



AT 
A2 



\BP J r \BP ) \BP 
If we contract the {xy} components of the equation of motion with Vy 

2 



B 

'a 



A' 



Q 



Q_ 

B^ 



AG A' 
3F \X 



^BP J ' \B^ 
On contracting with a unit vector orthogonal to Vy we obtain 

A' \ (A' 







S2 



IB 
+ 



BP 



BP 



52 ^ \By 



(198) 
(199) 

(200) 

(201) 

(202) 
(203) 

. (204) 

. (205) 
we obtain 
(206) 

(207) 
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12.1 Asymptotic behaviour 

We assume that, at large r, the leading asymptotic behaviour of the solution to 
these equations has the form 

^-^ + y + ^ + Oi^^) , (208) 

+ ^ + ^ + ' (209) 

and 

B^l + ^ + ^ + Oih . (210) 

On evaluating the two leading order contributions to each of the equations we 
find the results 

+ Si = . (211) 
This result holds in the conventional theory. 

A2 = ^A? + ^iAi , (212) 
A2 = \aI , (213) 

If we compare this asymptotic behaviour with the weak field approximation 
then we see that 

Ai^-GM , and A^^—FM , (215) 

where M is the mass of the object at the centre of the spherically symmetric 
system. The imphcation of these results for the barred (matter) metric is easily 
checked. We find 

^.v^i_2££ + P£i±|i£)M! , (216) 



and 



B^^l + Hl^^^ZiM . (217) 
r 



The standard post-Newtonian parametrisation is 



A2A2^1_^ + ?^^|^ , (218) 



and 

B^M^l + Mi^il , (219) 



24 



where for standard General Relativity we have /? = 7 = 1 . In the theory we 
have outlined we find 



/3 = 1 + |^-1 + |^ + 0(F^) , (220) 

and 



7 = 1-^ ■ (221) 

The time-delay results yield 7 = 1.00 ± .002 pT] , Q which places a strong re- 
striction on f /Gn- An even stronger restriction results from the very accurate 
VLBI data we have 7 = 1.000 ± .00032 It follows immediately that 

F < 2 X IO'^Gtv with the consequence /5 — 1 < 3 x 10"'^ . The appropriate 
conclusion is that on the scale oi G oy Gn-, F is very small. This however does 
not prevent it having an effect in the very early universe as indicated in the 
discussion of the FRW model. Nevertheless these experiments could be taken 
to rule out the model as we have presented it here. Variations of the model in 
which extra terms are introduced into the linking Lagrangian that give a mass to 
the new degrees of freedom will lead to them having a finite range. This would 
leave the long distance asymptotic behaviour of the model coincident with that of 
standard General Relativity at least in the weak field limit. It would still permit 
modifications of the standard theory at short distances. In particular it would 
allow the existence of very transparent modifications of the effective metric which 
could be interpreted as a novel form of matter. The potential of such a model 
to explain the presence of dark matter in the universe is intriguing and requires 
further study. 



13 Stability 

Because of the indefinite character of the local Lorentz metric it is not obvious 
that the degrees of freedom in the linking field M all contribute positively to the 
energy, in other words that they correspond to a stable theory. While we cannot 
at present give a final answer to this question we find it illuminating to consider 
the issue in the case of two dimensions. Here the matrix M has the form 




(222) 



with ad — he = 1 . The quantity / = TtMM^ is invariant under separate left 
and right Lorentz transformations. We have 

f = a^ + -b^ -c^ . (223) 

If we set 

M = LSL' , (224) 
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where L and L' are independent Lorentz transformations then of course 

/ = TrMM^ = Tr^^"^ . (225) 
If we choose the special form for S 

then 

/ = A^ + A-2> 2 . (227) 

If we choose 

y sm.9 cos 9 J ^ ' 

then 

/ = 2cos^^-2sin^^ = 2cos2^ . (229) 

In this case 

- 2 < / < 2 . (230) 

third possibility is 




(231) 



We then find 

/ = - (A^ + A~2) < -2 . (232) 

The three ranges together cover all possibilities for the value of / . Wc therefore 
can conclude that the matrix M will be gauge equivalent to the appropriate 
special form S provided that form is chosen according to the value of / . Indeed 
it is possible to check explicitly that given / the real matrices L, L' and S can 
be computed provided S is chosen in the appropriate manner. 

It is therefore clear that the object that results from equivalencing matrices 
M G SL{2,K) by simultneous left and right Lorentz transformations, although 
made up of pieces of manifold, is not itself a manifold. With the sign conventions 
we have adopted for the linking action, oscillations for which the field strengths 
maintain / > 2 will yield positive energies. The remaining two cases yield nega- 
tive and then positive energy oscillations. 

Similar yet more complicated possibilities exist for the four dimensional case. 
The completely positive mode case corresponds to matrices M that are Lorentz 
gauge equivalent to the special form S where 



/ Ao 











\ 





Ai 
















A2 






V 








As 


/ 



(233) 
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where A0A1A2A3 = 1 . This was the assumption we made imphcitly in the special 
cases we worked out above. It is not known whether solutions of the equations of 
motion remain in the purely in the positive sector. It is not known whether it is 
important for this to be the case. In the sense we have elucidated here, therefore, 
the stability of the theory remains undecided. However the examples discussed 
in detail show that possible solutions with stable characteristics exist. 



14 Conclusions 

We have shown that it is possible to formulate a purely geometrical theory with 
a variable speed of light in a covariant fashion. The theory gives a meaning to 
variations of the speed of light by identifying two vierbein bundles one associated 
with gravity and one associated with matter. Matter moves according a conven- 
tional Lagrangian in its own gravitational back ground. Its energy momentum 
tensor is conserved in the usual way. In this way it similar to other bi-metric the- 
10| , |Tl|, Gravity is also formulated in a conventional way in terms of 



its own vierbein bundle. The extra dynamics is supplied by a linking action that 
treats the matrix relating the two bundles as a dynamical variable. The form of 
the action is based on a sigma-model construction and requires the introduction 
of a new coupling constant with the same dimensions as the Newtonian constant. 
The action is not unique and can be modified by the addition of terms that ef- 
fectively give a mass to the new degrees of freedom at the expense of introducing 
new parameters. 

In this paper we have examined the the minimal version of the theory without 
extra mass inducing terms. As a result we find the asymptotic properties of the 
theory are modified relative to the standard Einstein theory. The parameter 7 
of the post-Newtonian approximation is not equal to unity. This is inconsistent 



with the time-delay and VLBI measurements |2J]. Since our theory 

reduces to the standard one when the new coupling constant, F, is set to zero 
we can of course achieve agreement with experiment by reducing this coupling 
to a sufficiently low value. A very low value for F does not by itself preclude 
strong effects in the very early universe. However the most natural conclusion 
is that the a departure from standard GR along the lines suggested is unlikely. 
It is possible that with the introduction of mass terms in the linking action the 
resulting modification of the asymptotic properties will leave the theory with 
essentially standard properties at long distances and with modifications only at 
short distances. We have not studied this point in this paper. 

We have also shown that a (fiat) FRW universe can be found in which it 
is possible to see the speed of light behaving in the way envisaged by previous 
authors. This would indeed permit the uniformisation of cosmic temperature 
associated with the cosmic microwave background. However we have not been 
able to analyse the model sufficiently closely to establish precisely what may 
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happen in the very early, or indeed late, universe. The results must therefore 
be treated as suggestive rather than definitive. Nevertheless the compatibility 
of the variable speed of light idea with covariance and energy and momentum 
conservation for matter has been established within a completely geometrical 
theory. 

Finally we must draw attention to the fact that the stability of the theory 
against the production of negative energy densities associated with the new de- 
grees of freedom has not been satisfactorily established. It would be interesting 
to check that matter oscillations do not give rise to negative energy radiation. 
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